Quantum states and linear response in dc and electromagnetic fields for charge current 
and spin polarization of electrons at Bi/Si interface with giant spin-orbit coupling 
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An expansion of the nearly free-electron model constructed by Frantzeskakis, Pons and Grioni 
[Phys. Rev. B 82, 085440 (2010)] describing quantum states at Bi/Si(lll) interface with giant spin- 
orbit coupling is developed and applied for the band structure and spin polarization calculation, 
as well as for the linear response analysis for charge current and induced spin caused by dc field 
and by electromagnetic radiation. It is found that the large spin-orbit coupling in this system may 
allow resolving the spin-dependent properties even at room temperature and at realistic collision 
rate. The geometry of the atomic lattice combined with spin-orbit coupling leads to an anisotropic 
response both for current and spin components related to the orientation of the external field. The 
in-plane dc electric field produces only the in-plane components of spin in the sample while both the 
in-plane and out-of-plane spin components can be excited by normally propagating electromagnetic 
wave with different polarizations. 
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I. INTRODUCTION 



It is well-known that the knowledge of a material with 
high values of spin-orbit (SO) coupling parameters is 
a goal for many theoretical, experimental, and device 
research groups in condensed matter physics and spin- 
tronics due to its fascinating spin-related properties and 
possible applications in the information processing and 
storage. Among the candidates which attract consid- 
erable attention throughout the last decade there are 
Bi/Si(lll) surface alloys which band structure has been 
experimentally studied for several years and recently 
modeled in a paper^ by Frantzeskakis, Pons and Gri- 
oni. This material, in line with the other examples of 
" metal-on-semiconductor" systems with large SO cou- 
pling, has been a subject of many both experimental 
and theoretical papers since it seems very promising to 
make use of a material combining the large SO split- 
ting of Bi and conventional semiconductor technology 
of Si which is one of the main goals of spintronicsi^^— 
Here we shall mention only some of the numerous re- 
sults of research on the Bi-covered Si surface proper- 
ties with various crystal orientation of Si substrate. In 
particular, the scanning tunneling microscopy has been 
used to determine the surface structure of Bi/Si some 
18 years ago,^ and the analysis of the atomic geome- 
try and electronic structure continued further^'^ focus- 
ing throughout the recent years mainly on the atomic 
surface geometry and spin-resolved band structure re- 
construction where the methods of angle-resolved pho- 
toemission spectroscopy have been applied"'^'^"^^. Other 
methods included the low-energy electron diffraction and 
atomic force microscopyiSiii, and besides pure Si, the 
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Si-Ge superlattices have been used as a substrate for 
Bi coverage^^ and later the lateral Ge-Si nanostructure 
prepared on the Si/Bi surface have been studied by the 
scanning tunneling microscopyJ^ Also, for the Bi/Si sys- 
tem there were studies of energetic stability and equi- 
librium geometry^ and the possibility of designing iron 
silicide wires along Bi nanolines on the hydrogenated Si 
surfacCj^i^ and of the thermal response upon the fem- 
tosecond laser excitation. It is well-known that Bi is 
a material with very big SO splitting, and thus it at- 
tracts a steady interest in its potential applications in 
spintronics where various schemes of combining it with 
semiconductors are suggested, one of the most recent be- 
ing an investigation of BiTel bulk material where the SO 
splitting reaches a very high value of 0.4 eVi^ 

From the list of papers mentioned above it is evident 
that the geometric properties of atom arrangement and 
the resulting band structure have already been studied 
for Bi/Si systems by many experimental and theoreti- 
cal groups. However, much lower attention have been 
given so far to the prediction and observation of different 
effects caused by the electron system response to an ex- 
ternal excitation, including such basic properties as the 
charge current and spin polarization in the dc field which 
are often are considered as the starting point of the re- 
sponse calculations, especially for systems with impor- 
tant role of SO coupling.— "-^X Besides the response to the 
dc electric field, the optical properties of SO-split band 
spectrum always attracted significant attention starting 
from the conventional semiconductor structures with big 
SO coupling. ^^"'^^ In our previous papers we have ob- 
served an important role of SO coupling in conventional 
InGaAs-based semiconductor superlattices on the energy 
band formation'^'^ which directly affected both charge and 
spin response for the excitation by the electromagnetic 
radiation^^i^ and by the dc electric field. '^^ It is known 
that the spin polarization configurations in semiconduc- 
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tors may have a rather long relaxation time^i^iiS^ which 
makes them as important as the conventional charge cur- 
rent setups for their applications in the nanoelectronics 
and spintronics. 

While there is no doubt that the electron properties 
in Bi-covered Si interface differ from the ones in conven- 
tional semiconductor structures with strong SO coupling, 
the questions regarding their SO-dependent response to 
the external dc and electromagnetic fields remain to be 
very important since we are still in the beginning of our 
way towards understanding and utilizing such novel ma- 
terials with giant SO coupling as Bi/Si. 

The goal of the present paper is to apply a modified 
and expanded version of simple but adequate nearly free 
electron model^ for the band structure of electrons in 
Bi covering the Si(lll) surface which allows us to cal- 
culate various physical characteristics of this material in 
the linear response regime, including the charge current 
and spin polarization caused by dc electric field with dif- 
ferent orientations, and also to obtain the response of 
non-equilibrium spin polarization excited by the electro- 
magnetic field with various polarizations. Indeed, such 
physical quantities can be among the first ones measured 
in the nearest experiments on Bi/Si, and hence it is of 
interest to calculate them beforehand both qualitatively 
and, when possible, quantitatively. Since we do not know 
exactly for the present moment many material parame- 
ters of the electron system in Bi-covered Si interface, in- 
cluding such parameters as the electron surface concen- 
tration and mobility, the dielectric tensor, the relaxation 
rates for charge and spin, etc, we sometimes cannot yet 
calculate the effects in the absolute measurable units and 
use the standard label "arbitrary units" instead. Still, 
we believe that the comparison of the output results for 
the same physical parameter calculated at different con- 
ditions always has a value since it allows to predict the 
relative significance of them when the conditions are var- 
ied. We find both common and distinct features of the 
charge and spin system response in Bi/Si compared to 
well-known GaAs or InGaAs semiconductor structures. 
Thus, we believe that our findings can be a good starting 
point for further theoretical and experimental studies of 
charge and spin response in Bi/Si system at various ex- 
ternal fields. 

The paper is organized as follows: in Sec. II we in- 
troduce the expanded version of the nearly free electron 
Hamiltonian for Bi/Si and discuss the band structure and 
spin polarization in the energy bands, in Sec. Ill we solve 
the kinetic equation and calculate the charge current and 
spin polarization in the dc electric field, in Sec. IV we 
study the excited spin polarization in the framework of 
the linear response theory for different polarizations and 
frequencies of the incident electromagnetic wave, and we 
give our conclusions in Sec.V. 



II. HAMILTONIAN AND QUANTUM STATES 

Recent proposal of several theoretical models for band 
structure calculation of Bi/Si electron surface states^ pro- 
vided a variety of choices for studies of the correspond- 
ing Hamiltonian and quantum states. Here we shall start 
with the simplest nearly free electron (NFE) model which 
was proposed initially for the description of the band 
structure mainly in the vicinity of the M point at the 
surface Brillouin zone (SBZ) of Bi/Si having a hexago- 
nal shape shown by dashed contour in FigHJ We shall 
only briefly describe it here since the detailed derivation 
and discussion is available in the original paperii The 
choice of the reciprocal lattice vectors initially has been 
restricted to three vectors Gi, G2 and G3 connecting 

four equivalent Gamma-points T^'^\ T^^\ F^^^ and T^^\ 
In the framework of the NFE approach for each Gamma- 
point the standard 2x2 Rashba Hamiltonian of a free 

electron in the basis has been written with the center 

— (") 

of the quasimomentum at the corresponding F point. 
As a result, an 8 x 8 matrix is derived giving the en- 
ergy bands and two-component eigenvectors (the Rashba 
spinors) describing the spin polarization in the reciprocal 
space. 

We are going to use an expanded version of this model 
by including the remaining reciprocal vectors G4, G5 
and Gg into our basis of nearest neighbor sites connect- 
ing the center Gamma-point F^''^ with all surrounding 

points f'"^'', . . . , f''^', as it is shown by hollow vectors 
Gi, . . . , Gg in FiglH where several hexagonal SBZ-s are 
shown by solid contours, thus creating a 14 x 14 Hamil- 
tonian matrix. We assume the previously determined^ 
values of geometrical parameters VM — 0.54 A and 
VK = 0.62 A . Such an expansion allows us to treat 
a much wider area of the SBZ compared to the region 
near the M-pointi and to keep the symmetry of the non- 
trivial hexagonal Bi/Si(lll) trimer structure with one 
monolayer of Bi atoms. ^—"i^ii^ 

Our Hamiltonian may thus be described via its matrix 
elements in the following form: 

Hnn' — H^Ck + G„)(5„„' 4- Vnn', (1) 

and the electron spinor wavefunction is constructed as 
^kW = ^ankV-nkW (2) 

n 

where the conventional form of Rashba Hamiltonian is 
used, 

fi - ( ^^fc^/2TO aR{ky + ikx)\ 
' ~ [ aniky - ik,) r^e/2m ) ' 

and the matrix elements of the periodic potential cou- 
pling the free electron states are 
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FIG. 1: Surface Brillouin zones (SBZ) of Bi/Si with a hexago- 
nal shape shown by sohd hnes and the reciprocal lattice (hol- 
low) vectors Gi , . . . , Ge connecting the equivalent Gamma- 
points r'°'', . . . , r''^' of the nearest neighbor approximation 
for the nearly free electron model. The spin-split Rashba 
paraboloids are centered in each of the Gamma-points. 



Vnn' = (V-n I ^ Vb exp(iG„r) I Vr,')- (4) 
m 

The basis functions V'nk(r) in Eq.(I2|) are the well- 
known Rashba spinors ■(/'nk ~ V'k+Gn where 

= _tg«Arg(fe^-ife,) ^ > (5) 

and the (±) sign corresponds to two eigenvalues for 
Rashba energy spectrum E{k) = ?i^fc^/2m ± ank. We 
continue to adopt here the known values of material pa- 
rameters for Eqs (H])-© and constants used during the 
initial construction^ of the NFE model for Bi/Si, namely, 
we put m — 0.8mo, — 1.1 eV • A and Vq = 0.3 eV. 

After diagonalization of the Hamiltonian ([T]) we ar- 
rive to the energy band spectrum E = Es{kx, ky) where 
s — 1,2,... labels the energy bands of the electrons in the 
Bi/Si(lll) system. The three-dimensional plot of the en- 
ergy band structure is presented in Figl2]for the four low- 
est bands labeled from 1 to 4. These lowest bands seem 
to be of the primary importance for the electron response 
analysis since the Fermi level is reportedly located!, in the 
middle of them at Ep = 1.6 eV, i.e. between the band 
No. 2 and band No. 3, as it can be seen in Figl2l One of the 
most important features of the spectrum in Figl2] stem- 
ming from the lattice geometry is the hexagonal symme- 
try of the energy bands in k-space which implies, among 
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FIG. 2: (Color online) Energy band structure of the electrons 
on Bi-covered Si(lll) surface for the four lowest bands labeled 
from 1 to 4. The Fermi level is located at Ef ~ 1.6 eV 
between the band No. 2 and band No. 3 where a small global 
energy gap of around 0.1 eV is formed. The cross-sections of 
our 3D band plot shown here accurately repeat the 2D plots 
for the energy dispersion lines along various directions in the 
SBZ which were studied earlier in the framework of the NFE 
model.— 

other things, the absence of the -O- ky symmetry, lead- 
ing to rich properties of the spin response phenomena as 
we shall see below. It should be mentioned that the cross- 
sections of our 3D band plot shown here accurately repeat 
the 2D plots for the energy dispersion lines along vari- 
ous directions in the SBZ which were studied earlier in 
the framework of the NFE modelii We would like to add 
here just some new quantitative data: the cross-sections 
of the energy band surfaces along the T — M direction 
reported previously^ might created an impression of a 
large energy gap in the whole spectrum between bands 2 
and 3 while the complete 3D presentation of these bands 
in Fig l2] indicates that the global structure of the energy 
bands in the whole 2D SBZ leaves this gap opened but 
with a much smaller width of around 0.1 eV. Of course, 
the precise values of energy gaps may vary from model 
to model and can be specified more precisely during the 
future experimental and theoretical analysis. 

Another important characteristic of quantum states in 
a system with significant SO coupling is the spin polar- 
ization of the eigenstates ipk in the Brillouin zone which 
can be defined as the vector field in the reciprocal space 
with the components (m = x, y, z) 

S'„i(k) = (l/'k I CTm I Ipk)- (6) 

As usual for the Hamiltonian with a pure Rashba SO 
coupling term, the out-of-plane component Sz of the spin 



FIG. 3; Two-dimensional spin polarization distributions 
(51,(14), Sy(k)) for the lowest energy bands No.l (a) and No. 2 
(b) from the band spectrum shown in Fig[2l with the hexag- 
onal SBZ marked by a solid contour. The initial Rashba 
counter clock-wise and clock-wise patterns are present in a 
wide area surrounding the SBZ center, but more complicated 
vector field structure near the SBZ edge. 



field vanishes. The remaining components form a 2D spin 
polarization distribution in the SBZ which forms a spe- 
cific vector field picture for each of the energy bands. In 
FigEland Fig|4]we show the 2D spin polarization distri- 
butions {Sx(k), SyCk)) for the lowest energy bands No.l 
(Fig Hi), No.2 (Figll)), No.3 (FigUi) and No.4 (Figgb) 
from the band spectrum shown in Figl2l with the hexag- 
onal SBZ marked by a solid contour. As for the spins 
in two lowest subbands shown in FiglSj one can see here 
that the initial Rashba counter clock-wise and clock-wise 



FIG. 4: Spin polarization distributions {SxO<.), Sy{k)) for the 
higher energy bands No.3 (a) and No.4 (b) from the band 
spectrum shown in Fig|2] with the hexagonal SBZ marked by 
a solid contour. The spin polarization in band No.3 (a) and 
especially in band No.4 (b) demonstrates more new properties 
compared to the free Rashba states, including both the shape 
of the spin vector field which captures now more features of 
the hexagonal geometry of the SBZ, and the arising of new lo- 
cal vortices at various points of symmetry of the SBZ, mainly 
near its corners. 



patterns of spin directions are present in a rather wide 
area surrounding the SBZ center, but more complicated 
vector field structure arises near the SBZ edge. The spin 
polarization in higher band No.3 and especially in band 
No.4 shown in Fig|4] demonstrates more new properties 
compared to the free Rashba states, including both the 
shape of the spin vector field which captures now more 
features of the hexagonal geometry of the SBZ, and the 
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arising of new local vortices at various points of symme- 
try of the SBZ, mainly near its corners. 

It can be concluded from the analysis of the spin polar- 
ization in the energy bands of Bi/Si system that certain 
properties of the initial basis of Rashba states remain vis- 
ible. However, the change of the space symmetry to the 
hexagonal type without the element of axial symmetry 
and without the x ^ y symmetry may probably lead to 
both common and distinct features in the current and 
spin response to the application of various external fields 
compared to the well-known properties of 2DEG with 
Rashba SO coupling. This assumption will be confirmed 
and illustrated below. 



III. CHARGE CURRENT AND SPIN 
POLARIZATION RESPONSE FOR DC FIELD 

It is known that the response of a two-dimensional elec- 
tron gas with SO coupling to a constant electric field 
may be accompanied not only by the charge current but 
also by the spin polarizationi^i22i2& The most significant 
properties of such response for a pure Rashba SO term (|3]) 
in the Hamiltonian is the arising of the in-plane trans- 
verse polarization, i.e. the Sy(^^) spin component when 
the electric field is applied along the x{y) direction while 
the out-of-plane Sz component is absent in case of the 
accurately included relaxation processes which is some- 
times related also to the absence of the spin Hall effect for 
a k-linear Rashba model in the presence of the disorder.— 
So, it is natural to start the analysis of the electron sys- 
tem response with the calculation of the charge current 
and field-induced spin. 

We shall start with the calculation of the non- 
equilibrium electric field-affected distribution function 
/m(k, Ei) in the m-th energy band. If the system is sub- 
jected to a constant and uniform electric field Ei parallel 
to the i-th axis, then in the collision frequency approxi- 
mation the kinetic equation for /m(k, Ei) can be written 
as^ 

- -^[/™(k,i?.) - /„.(k)], (7) 

where i' is the collision rate and /m(k) = 1/(1 + 
exp[(£'„j(k) — /i)/fcBT]) is the Fermi equilibrium distribu- 
tion function in the m-th band. Since the Bi/Si energy 
spectrum is characterized by a very large SO splitting 
and the band widths of the order of 1 eV, it may be a 
promising candidate for spin-dependent phenomena visi- 
ble at room temperature. Thus, in the following we shall 
assume that T = 293 K and consider a value of = 10^^ 
s~^. As we have said before, the estimate for the col- 
lision rate as well for many other material parameters 
for the Bi/Si system is presently based on the assump- 
tions rather than of the solid experimental facts since 
we are still in the beginning of the investigation for this 
new material. Still, we believe that our qualitative and 



sometimes quantitative results can be useful for predict- 
ing some novel properties of the electron and spin system 
response. 

The mean charge current density ji{Ei) measured for 
2D system in units of current divided by the unit of 
transverse system size and the mean spin polarization 
Sj, j — X, y, z, can be found as^ 

ji{Ei) = en^(i/'mk I I i/'mk)/m(k, E'l), (8) 

S-AEi) = ^(^Z-^k I <Jj I V^mk)/m(k, Ei), (9) 
m,k 

where n is the surface concentration of the electrons 
on the Bi-covered Si surface, Ui {i = x, y, z) are the Pauli 
matrices, and the velocity operator Vi = dH/dki includes 
the SO part proportional to the Rashba parameter an 
and acts on the spinor wavefunctions ([2]) via the matrices 

( -insiy/m a/n \ 
"^v- \ a/n -invy/m ) ■ 

In order to take a closer look for onto the expectations 
of the charge current density, we shall calculate jx{Ex) 
and jy{Ey) as well as Si{Ex) and Si{Ey) {i — x,y,z) 
by applying Eq.® for the electron surface concentra- 
tion n — 10^^ cm^^ which may be reasonable since the 
surface density of atoms on Bi-covered Si(lll) surface 
is estimated on the order of 10^^ cm~^ according to the 
experimental observations;^^ 

The results for the charge current density (|8]) and spin 
polarization ^ are shown in FiglSJa) and Figl5l[b) re- 
spectively. One may see a conventional linear dependence 
of the charge current on the applied electric field through- 
out the whole range of fields up to 2 kV/cm, and almost 
linear dependence for the significant non-zero induced 
spin components Sy{Ex) and Sx{Ey). Other in-plane 
components Sx{Ex) and Sy{Ey) marked by arrow are 
also present in Fig|S] but their magnitude is much lower 
compared to Sy{Ex) and Sx{Ey), and the out-of plane 
Sz component is negligibly small. It is evident that the 
lattice asymmetry with respect to the x y interchange 
transform has lead to a slight but distinct asymmetry in 
the current amplitude of around 12%, and the dominat- 
ing Sy{Ex) and Sx{Ey) induced spin components demon- 
strate the well-known transverse in-plane character of the 
induced spin for linear Rashba SO coupling. It should be 
noted that the local probe measurements of induced spin 
polarizations (or magnetization) may detect the non-zero 
static and dynamic local magnetization^S. in the spot un- 
der the probe even in case of total zero mean spin value 
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FIG. 5: (a) Charge current density ^ and (b) spin polar- 
ization ((9)1 induced by the dc electric field applied along x 
(solid curves) and y (dashed curves). A conventional linear 
dependence of the charge current on the applied electric field 
throughout the whole range of fields and almost linear depen- 
dence for the significant non-zero induced spin components 
Sy{Ex) and Sx{Ey) are visible. Other in-plane components 
Sx{Ex) and Sy{Ey) marked by arrow also present but their 
magnitude is much lower, and the out-of plane Sz component 
is negligibly small. 



IV. SPIN POLARIZATION EXCITED BY 
ELECTROMAGNETIC FIELD 



The response of the spin system in materials with sig- 
nificant SO coupling on the application of an external 
electromagnetic radiation is among the most important 
and straightforwardly obtained characteristics since the 
optical manipulation of spins is one of the main goals of 
spintronics in general, and the linear response theory for 
the electromagnetic radiation effects is well-established 
and easily applied. It was found in various papers that 
the spins with different polarizations can be excited, de- 
pending on the symmetry of the electron Hamiltonian, 
the type and strength of the SO terms, and on the po- 
larization of the incident radiationi^i^i^^"— i^ii^ As in the 
previous Sec, we shall calculate the response functions 
for the room temperature and for a realistic collision 
broadening since the relatively large scale of energy bands 
and SO splitting in the Bi/Si electron system compared 
to the conventional GaAs, InGaAs or pure Si semicon- 
ductor structures can make Bi/Si being a promising can- 
didate for the observation and control of the predicted 
radiation-induced effects in the devices operating even at 
room temperatures, as we hope. 

The electromagnetic wave is considered to be propa- 
gating normally to the Bi/Si(lll) interface along the z 
axis, and is characterized by the polarization of the elec- 
tric field vector E = Eq exp(i(k-r— wt)) in the (xy) plane, 
Eo — {EQj:,Eoy,0). In the dipole approximation the in- 
teraction of the electromagnetic field with the electrons 
is described via the velocity operators (ITU|) . p^ which in- 
clude the SO part. We start with the calculation of the 
absorption coefficient 



The examples of such systems with zero total spin 
polarization but non-zero spin spatial distribution (spin 
textures) can be found among the models of semicon- 
ductor superlattices with SO coupling both with2^ and 
withou t'^^'^^ external magnetic field, but their experimen- 
tal observation and device application are currently lim- 
ited by the probe and manipulation technology of the size 
of artificial superlattices and quantum wells rather then 
probing and utilizing the spatial magnetic configurations 
on the scale of individual atoms in the lattice. As to 
the predicted non-zero mean values of the induced spin 
such as those predicted here for Bi/Si, they are related to 
the whole sample and thus should be detectable. We be- 
lieve that the predictions of the charge current and spin 
polarization generation made in this Sec. can be useful 
in designing novel spintronic devices where the induced 
spin components are coupled in a well-defined manner to 
the direction of the applied electric field, and this effect 
survives at room temperature and finite collision rate. 



a{u^) = — — ^ I (e • v)„„,^ p X (12) 

xS{En'k - Enk ~ huj) ifnk - fn'k) (13) 

where e = (bx, By, 0) is the polarization vector for the 
incident wave, v = {vx,Vy,0) is given by the velocity 
operators (fTO |) . (|ll|) . /„k; and /n'k are Fermi equilibrium 
distribution functions, V is the sample volume, and the 
summation is taken over all energy bands n, n' and the 
SBZ points k. As we have already said, we don't know 
yet the exact values for many of the material param- 
eters for Bi/Si including the dielectric constant e, and 
thus we shall focus mainly on the dependence of ([TS]) on 
the incident wave frequency and will scale the absolute 
value of absorption in arbitrary units which can always 
be rescaled when the values material parameters will be 
clarified in future experiments. 

The second quantity which frequency dependence we 
shall present together with the absorption is the induced 
spin polarization Sm{^) which can be derived by apply- 
ing the Kubo linear response theory^^i^i^ and has the 
following form: 
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'(k) 



(14) 
(15) 



Here the interband matrix elements of the spin m-th 
component operator sl^^ik) = {ipn'k \ \ ipnk) as 
well as the matrix elements for the ^-th component of 
the velocity operator p^ . ([TT|) enter depending on the 
incident wave polarization and on the desirable output 
for the spin component, ?iw„„'(k) = i?„k — £'n'k, the 
parameter v is level broadening which we take as being 
equal to the collision rate introduced in the previous Sec, 
and Vc is the unit cell volume. As before, we shall assume 
that T = 293 K and ly = 10^^ which should provide 
us with realistic expectations regarding the absorption 
and induced spin dependence on the incoming photon 
energy. 

The results for the photon energy dependence of the 
absorption coefficient (|13p and the induced spin polariza- 
tion dTSl) arc shown in FiglS]through Figl5]for the incident 
wave linearly polarized along x (Fig|6]), along y (Fig 17]) 
directions, and for circular cr+ (FiglSJ and cr+ (FiglH]) po- 
larizations. The photon energy interval is chosen to cover 
the whole energy band range of the four lowest bands 
shown in Figl2] where the most effective transitions oc- 
cur between the states below and above the Fermi level. 
One can see that both the in-plane spin components Sx, 
Sy and the out-of plane component Sz can be excited on 
a comparable scale which is a distinguishable feature of 
the lower hexagonal symmetry combined with Rashba SO 
coupling compared with one-dimensiona l^^'^^i'^^ or two- 
dimensional square2S lattices with Rashba SO coupling. 

If one compares the results for x- and y-polarized radi- 
ation in Figiniand Figl?! it can be seen that, similar to the 
dc current properties discussed in the previous Sec, the 
X y lattice and energy band asymmetry in the hexago- 
nal geometry of the whole problem is reflected here in dif- 
ferent shape and amplitude for the absorption coefficients 
in Figini^a) and[71[a). Again, by looking onto the relative 
magnitude of the excited spin components in Fig|6l[b)-(d) 
and in Figl7Kb)-(d) one can see that the pure Rashba SO 
coupling is reflected in the dominating excited Sx compo- 
nent in Figini and correspondingly in the dominating Sy 
component in Fig|6j i.e., in the in-plane and transverse 
direction relative to the electric field vector of the inci- 
dent wave. As to the results for the circular cr^-polarized 
radiation presented in FigjS] and FiglHl on can see that 
the common features of the response to both a;-polarized 
and y-polarized incoming wave from Fig 15] and Fig |7] can 
be seen on the spin component figures since both Vx and 
Vy operators here enter the expressions (IT^ and (|T5]) for 
the response. The absorption coefficient is totally insen- 
sible to the direction of rotation for the incoming wave as 
it can be seen by comparing FiglS][a) and Figin][a). The 
shape of the photon energy dependence for the excited 



(^a) A Absorption (arb. units) 
0.08 - 




0.5 1.0 

A S (arb. units) 
0.2P ^ 



h(jd (eV) 

A 




-0.3 - 



(d) A ^ (arb. units) 



h(A) (eV) 




FIG. 6: (a) Absorption coefficient and (b) - (d) components 
of the induced spin density in the Bi/Si surface electron gas 
shown as a function of the incoming photon energy for the 
linear a;-polarized incident radiation propagating normally to 
the interface plane. The greatest absorption and most of the 
spin polarization peaks are induced in the photon energy near 
Twj ~ 0.5 eV corresponding to the transitions between the 
highest occupied band 2 and the lowest unoccupied band 3 
of the electron energy spectrum shown in Fog|5] Both the in- 
plane spin components Sx , Sy and the out-of plane component 
Sz can be excited. 
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FIG. 7: Photon energy dependence of (a) absorption coeffi- 
cient and (b) - (d) components of the induced spin density 
for the same parameters as in Fig[6] but for the y-polarized 
incident radiation. The x ^ y lattice and energy band asym- 
metry in the hexagonal geometry is reflected in different shape 
and amplitude for the absorption coefficients in this Fig. and 
in Fig|6l The Rashba term in SO coupling is reflected in the 
dominating excited Sx component in this Fig. and the domi- 
nating Sy component in FiglU 
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FIG. 8: Photon energy dependence of (a) absorption coeffi- 
cient and (b) - (d) components of the induced spin density 
for the same parameters as in Fig[6]but for the circular a^- 
polarized incident radiation. The common features of the 
response to both a;-polarized and y-polarized incoming wave 
from Fig|6]and Fig[7]can be seen on the spin component fig- 



9 



(a) A Absorption (arb. units) 




ha> (eV) 



^, 0.5 1.0 

WA 5, (arb. units) 



1.5 



2.0 2.5 



TkjO (eV) 




Cc)^^S^(arb. units) 



h(jO (eV) 




FIG. 9: Photon energy dependence of (a) absorption coef- 
ficient and (b) - (d) components of the induced spin den- 
sity for the same parameters as in Fig[6]but for the circular 
(T~ -polarized incident radiation. The absorption coefficient 
is unchanged from the cr^ case in Fig|8]^a) while the excited 
spin components demonstrate quantitative differences while 
maintaining the same type of shape for the photon energy 
dependence. 



spin components in Figl8l^b)-(d) and Figl9l^b)-(d) is dif- 
ferent for (T+ and a~ polarizations, but these differences 
have a quantitative rather than a qualitative character 
since the hexagonal symmetry of the lattice and the en- 
ergy bands does not make any of these two polarizations 
preferable from the point of view of the response quanti- 
ties (Uni) and (fT5l). 



In conclusion, the calculation and analysis of the ab- 
sorption and spin polarization response to the monochro- 
matic electromagnetic radiation with normal incidence 
and having different polarizations demonstrates that this 
radiation is most effectively absorbed in the photon en- 
ergy range of around 0.5 eV corresponding to the pho- 
ton wavelength A = 2.47 /J,m where both in-plane and 
out-of-plane spin components can be excited at realistic 
temperature and collision broadening on a comparable 
scale with relative amplitudes depending on the precise 
value of frequency and the polarization of the incident 
radiation. These properties can be useful for designing 
new optical and spintronic devices coupling the electron 
spin with light and operating at room temperature. 



CONCLUSIONS 



We have developed an expansion of the nearly free- 
electron modeli describing the energy bands and spin 
polarization for the electron states at Bi/Si(lll) inter- 
face with giant spin-orbit coupling, and applied it for the 
linear response analysis for charge current and induced 
spin caused by dc field and by electromagnetic radiation. 
It was found that the large spin-orbit coupling in this sys- 
tem may allow resolving the spin-dependent properties 
even at room temperature and at realistic collision rate. 
The geometry of the atomic lattice combined with spin- 
orbit coupling leads to an anisotropic response both for 
current and spin components related to the orientation of 
the external field. The in-plane dc electric field produces 
only the in-plane components of spin in the sample while 
both the in-plane and out-of-plane spin components can 
be excited by normally propagating electromagnetic wave 
with different polarizations. The qualitative predictions 
of the charge and spin response in a novel and promising 
Bi/Si system may be useful for the forthcoming detailed 
theoretical and experimental studies which may lead to 
the development of principally new electronic, optical 
and spintronic devices operating at room temperature. 
Further theoretical and especially experimental studies 
of this promising system with big SO coupling allow- 
ing the survivability of the spin-related effects at room 
temperature are expected bringing us new and fascinat- 
ing phenomena with both fundamental, experimental and 
device-related results. 
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